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Cell dynamics simulation is used to investigate the phase behavior of block copolymer/homopolymer mixture 
subjected to a steady shear flow. Phase transitions occur from transverse to parallel and then to perpendicular 
lamellar structure with an increase of shear rate and this is the result of interaction between the shear flow and 
the concentration fluctuation. Rheological properties, such as normal stress differences and shear viscosity, are 
all closely related with the direction of the lamellae. Furthermore, we specifically explore the phase behavior 
and the order parameter under weak and strong shear of two different initial states, and realize the importance 
of the thermal history. It is necessary to apply the shear field at the appropriate time if we want to get what 
we want. These results provide an easy method to create ordered, defect-free materials in experiment and 
engineering technology through imposing shear flow. 
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1. Introduction 

One of the key targets in modern material science is preparation of highly ordered and controllable 
nano-structures. Spontaneously formed structures usually do not exhibit long-range order, and often con¬ 
tain a large number of defects, such as dislocations and disclinations. Thus, it is an important prerequi¬ 
site that the emerging structures should be relatively defect-free. The long-range order can be induced by 
external guidance such as the shear field due to the inherent softness of polymers. The self-assembly phe¬ 
nomena of polymers under shear field analogous to those in a range of synthetic and biological systems 
can provide routes to the creation of novel materials. Furthermore, it is also easy, simple and operable to 
apply the shear field. Thus, this is a versatile means to obtain the long-range order and then create mi¬ 
crostructures with potential applications in biomaterials, optics, and microelectronics. The use of shear 
has proven to be an excellent method for achieving long-range order during the past decades, and the 
phase behavior and rheological properties of polymeric materials subjected to shear flow also has made 
a remarkable progress in both experiments and simulations fltilTl. Earlier, the study on block copolymer 
concentrated on experimental investigations, and most of them observed various alignments in a lamel¬ 
lar diblock copolymer by using TEM and SAXS, polystyrene-polyisoprene being the most representative. 
Parallel and perpendicular orientations were observed by Winey et al. in 1993 (12l and Patel et al. (13) 
in 1993, respectively; transverse orientation was found by Gupta et al. (14) and Zhang (l3, and it was 
defined as the one having a lamellae oriented perpendicular to both the shearing surfaces and the shear¬ 
ing direction. Then, Pinheiro and Winey observed mixed parallel-perpendicular morphologies in diblock 
copolymer systems at intermediate temperatures (16). Afterwards, researchers theoretically predicted a 
phase transition of the lamellar I17tl20l . ring (21), hexagonal cylinders (22), diblock copolymer subjected 
to shear flow further. These predictions were proved by researchers using various numerical simulation 
methods, such as ceU dynamics simulation (CDS) 12311291 . nonequilibrium molecular dynamics simulation 
(NEMD) I291132I . dissipative particle dynamics method (DPD) I331t35l . self-consistent field theory (SCFT) 
and lattice Boltzmann (LB) method (36), Brownian dynamics (BD) (37), mean-field approach (MFA) 1^ . 
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density functional theory (DFT) (3^|40), Monte Carlo simulation (MC) (42, molecular dynamics simu¬ 
lation (MD) (13 etc., Luo et al. discovered that the parallel orientation is stable at low shear rate and 
the perpendicular orientation is stable at high shear rate in the hexagonal cylinder phase of asymmet¬ 
ric block copolymer, and explored the kinetic pathways based on the time-dependent Ginzburg-Landau 
approach (26). Later on, the phase transition of gyroid phase (27), and sphere phase (28] of an asym¬ 
metric block copolymer under shear was also investigated. As for the symmetrical block copolymer, Lisal 
and Brennan’s simulation results indicated that the perpendicular lamellar phase persists for all shear 
rates investigated, whereas the parallel lamellar phase Is only stable at low shear rates, and it becomes 
unstable at high shear rates (^. In addition, Guo studied the amphiphUic model system: Its kinetics of the 
shear-induced isotropic-to-lamellar transition (^, and the paraUel-to-perpendicular orientation transi¬ 
tion in the amphiphilic lamellar phase under shear flow (30||33. Fraser et al. investigated how shear flow 
affects the orientation of lamellar structures. They found out that for any finite shear rate, the entropy 
production and the average energy are the smallest if the lamellae’s normal is perpendicular to the shear 
and the shear gradient direction, and this finding strongly suggests that the orientation corresponds to 
the global minimum for the bulk 02]. For all these cases, the shear flow plays an important role as a 
means for aligning the microscopic domains. 

Even though the phase behavior of diblock copolymer systems under shear flow Is relatively well 
understood, only a limited amount of information is available on the response of a complex polymer 
to shear. However, during the polymer processing, it is always doped with other substances in a block 
copolymer with the development of Industrial synthesis technology, such as homopolymer, nanoparti¬ 
cle. The rheological behavior of binary mixtures of a polystyrene-polyisoprene-polystyrene (SIS) copoly¬ 
mers and homopolymer polystyrenes (PS) were studied by Baek and Han (43). It is showed that the 
rheological behavior of mixtures is closely related to their morphological state, which in turn depends, 
among many other factors, on the ratio of the blends, molecular weight of the added homopolymer, and 
temperature. Further experiments also Investigated the shear-induced phase behavior of block copoly¬ 
mer/homopolymer (21 and block copolymer/epoxy (45) blends. Later on, some researchers focused their 
attention on the phase behavior of block copolymer/nanoparticle 14611511 composites under shear flow. 
Hence, He et al. systematically investigated the shear-induced reorientations and phase transitions of 
symmetric diblock copolymer/nanorod nanocomposites subjected to a steady shear flow l50l and oscilla¬ 
tory shear flow (2) via DPD. 

Consequently, it is really of great importance to explore the mysteries of block copolymer/homopolymer 
under external fields. Thus, the rest of this paper is devoted to the investigation of the phase behavior and 
rheological properties in a block copolymer/homopolymer mixture by using the cell dynamics simulation 
of time dependent Ginzbrug-Landau theory proposed by Oono and co-workers I52ii56l . It is expected that 
this will provide some guidelines for experimentalists. The other parts of this paper are organized as 
follows: section 2 is devoted to the description of the model and simulation method; section 3 is the nu¬ 
merical results and discussions; and finally, section 4 gives a brief conclusion of this work. 


2. l\/lodels and simulation methods 

We consider a block copolymer/homopolymer mixture rapidly quenched into the spinodal region, 
then subjected to a steady shear flow. Each copolymer chain is composed of monomers A and B with a 
short-range repulsive interaction between them. The interaction between C monomer in a homopolymer 
and B monomer in a copolymer is also assumed to be repulsive to each other. Hydrodynamic effects that 
prevail at the very late stage of phase separation in a polymer blend are neglected in the present model. 

Several parameters are defined to describe the system. We take (pB 0 > and (pco as the average 
volume fractions of monomers A, B, and C. In the case of symmetric block copolymers which is con¬ 
sidered here, polymerization indices of the A and B blocks are the same; this is ensured by - (pBo- 
An important quantity, the composition ratio of homopolymers and copolymers which can be changed, 
is defined by / = fco/ifAO + /bo) - In the process of phase separation, fluctuations are dominant, so 
we should investigate the local volume fractions of monomers A, B, and C. They are denoted, respec¬ 
tively, by (pAix, y, z], (psix, y, z], (pcix, y, z). Under the incompressibility condition, that is the total density 
(pA^x,y,z) -I- (pBix,y,z] + cj)clx,y,z] is constant, two of the local volume fractions are independent. Then, 
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we take (p[x, y, z) - (pA^x,y, z] - (ps{x, y, z), and i/r(x, y, z) - (pAix, y, z) + (pB(x,y, z) are independent vari¬ 
ables used to characterize the structure ordering. The order parameter (p{x, y, z) gives the local concentra¬ 
tion difference between the A and B monomers, the order parameter \f/{x, y, z] describes the segregation 
of the homopolymer and the copolymer. 

Here, we use a three-order-parameter model In (^. The free-energy functional of the system is given 
by 

F^Fl + Fs, ( 2 . 1 ) 

the long-range part Fl and the short-range part Fs are given by 

Fl= j JJ drdr'G[r,r')[(p(r)-(po]l(p(r')-(po], (2.2) 


and 


Fs 


= ^dxdy|y(Vi//-)^-i-y(V(/))^-i-w(i//,(/))], 


(2.3) 


respectively, where a, ci, cz are all positive constants. The long-range part is relatively simple, in which 
G{r,r') is the Green’s function defined by the equation -V^G(r,r') = 8[r- r'), while </)o is the spatial 
averages of (p. We should set (po-Oin the case of symmetric copolymers. As for the short-range part, the 
Cl and C 2 terms correspond to the surface tensions. The local interaction term LD(y/,(p] could be replaced 
by aj{ri,(p), where rj - y/ - y/c with y/c being the volume fraction at the critical point of the macrophase 
separation and being a constant determined by the parameters of the system. According to Ito et ah, 
we obtain 


a){ri,(p) ^ viir]) + vzicp) + biT](p- 



(2.4) 


where the functions viirj), vzicp) are assumed to be even with respect to the arguments, hi, bz are all 
positive constants. Equation {2A) prescribes a minimal model of the short-range part of the free energy 
in our system. In the symmetric case, hi = (-XAC + XBc')l2,bz- l/(i//’^ArA). It is now clear that the hi term 
mainly originates from the interaction between the A-B copolymer and the C monomer. If the repulsive 
interaction strength xbc between B and C is large enough, hi is positive, bz term represents the fact that 
the microphase separation should occur only in the copolymer-rich phase. 

In terms of the free energy functional in equations we obtain a set of two coupled equations 


dr) 

— = M„V 
dt 


.6F{i],(P) 

8rj 


-v-Vrjir, t), 


(2.5) 


d(p 


>8FIt],cP) 

8(p 


-v-V(p{r, t), 


( 2 . 6 ) 


where and Mip are transport coefficients, v is an external velocity field describing the shear flow 
profile 

v^fyex, (2.7) 


where f Is the shear rate and ex is the unit vector in the x direction. Numerical solutions of the above 
model system can be carried out in an LxX Ly x Lz three-dimensional cubic lattice, by using the cell 
dynamics simulation (CDS) approach, mainly because It is a fast method to simulate kinetic processes 
in phase separating systems of large sizes. The order parameters for each cell are rjln, t], y/{n, t), where 
n = {rix, riy, riz) is the lattice position and rix, riy, and riz are integers between 1 and L. The CDS equations 
corresponding to equations ^2.5M2.6) , in their space-time discretized form, are written as follows: 


rj(x,y,z, t+ At] 
(p[x,y,z, t+ At] 

where 


riix, y, z, t] + (((/^» - 1^] “ [v^x +l,y,z,t]-r][x-l,y,z,t]], ( 2 . 8 ) 

cp(x, y, z, t] + {{(1(1,)) - Iff,) - ^fy [(p[x +l,y,z,t]- (p{x - 1, y, z, r)], (2.9) 


Irj = -Di[{{T]))-r])-ArjtanhT] + r]+bicp-^bzT](p, (2.10) 

ly, = -Dz{{{(p))-(p)-A(i,tanh(p + (p+biT]-bzr](p, (2.11) 
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and 


6 ^ 3 ^ 1 

<<x» =—) x +—) x+ — 
80 ^ 80 80 


E 

NNN 


( 2 . 12 ) 


the subscripts N, NN, and NNN stand for the nearest-neighbor, the next-nearest neighbor, and the next- 
next-nearest neighbor cells, respectively To perform the numerical operations in a cell dynamics 

system, in general, the lattice size (Ax,Ay or Az) and the time step (Af) are all set to unity. 

We choose x-axis as the flow direction, y-axis as the velocity gradient direction and z-axis as the 
vorticity axis. A shear periodic boundary condition proposed by Ohta et al. 15911621 should be applied to 
X direction, and the periodic boundary conditions are applied in the y and z direction. With the shear 
strain y, this boundary condition is written as follows: 


T]{nx, riy, Hz, t)-T][nx + NxL + jWNyL, riy + NyL, riz + NzL] , (2.13) 

(p{nx, Uy, Hz, t)^(l)[nx + NxL + j{t)NyL, Uy + NyL, riz + NzL] , (2.14) 

where Nx, Ny, and Nz are arbitrary integers. 

For the three-Dimensional space, the first and the second normal stress differences Ni and A/ 2 , are 
defined as (3]|63l 

M = -|^[fc'-fcy]s(fc,r), (2.15) 

Nz = -J ^ [k^-kl] S[k, t), (2.16) 

and the excess viscosity defined as 16411661 

1 r dk 

—^kxkyS(k,t). (2.17) 

7 J [27T)-^ 


Our simulations are carried out on Lx x Ly x Lz - 64 x 64 x 64 cubic lattice, and the parameters are 
chosen to be Arj - 1.3, = 1.1, Di = 1.0, Dz - 0.5, hi = 0.10, bz - 0.02, and = Mfj, = 1, the initial distri¬ 

butions of (p and 77 are specified by random uniform distributions in the range [- 0 . 01 , 0 . 01 ], in accordance 
with the previous work 16711691 . In this paper, all parameters are scaled, so all of them are dimensionless 
(52). In the present simulation work, a - 0.02 and the formula is a = 12 / [AT^/b (1 - /b)] (/b denotes a block 
ratio) rTOl . 


3. Numerical results and discussion 

In order to investigate the phase behavior and rheological properties of the block copolymer and ho¬ 
mopolymer mixture under steady shear, CDS simulation is performed on the systems. The shear direction 
and the schematic of three different lameUar orientations are shown in figure and they will be men¬ 
tioned in the following discussion. Under zero-shear, we obtain abundant morphologies with different 
composition ratios of copolymer and homopolymer that are not displayed in this section. They conclude 
transverse, parallel, perpendicular lamellar structure and the concentric cylindrical structure that the 
block copolymer is surrounded by homopolymer when the homopolymer is in majority. Next, we want 
to know the structures corresponding to different shear rates relative to the transverse lamellae under 
zero-shear. Here, / = 55/45, all other conditions are the same except the shear rate. It means that the 
initial state at all shear rates includes a zero-shear and is disordered. The imposed shear rates 7 range 
from 0 to 0.001 and we show a part of them in figure]^ 

The ordered transverse lameUar structure mentioned above under zero-shear is shown in figure[^(a). 
With an increase of the shear rate, the chain stretches along the x-axis direction until it stretches com¬ 
pletely. Owing to a different intensity of the shear flows, we obtain different domain morphologies. The 
structure shows that the width of the domain is larger than the others corresponding to the shear rate 
7 = 0.000001 in figure[^(b), because the shear rate is too small to stretch entirely. We need to consider the 
orientations of the model after the chain stretches completely along x-axis at an appropriate shear rate. 
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Figure 1. Schematic representation of the geometry considered including the shear direction and three 
different lamellar orientations discussed in the text. 



(d) (e) (f) 


Figure 2. (Color online) Block copolymer/homopolymer under different shear rates, y = 0.0,0.000001, 
0.000005, 0.000006, 0.0006, and 0.001 in (a)-(f). Phase A is represented by the red regions, phase B by the 
yellow regions, phase C by the blue regions. The arrows in (a) show the direction of the x-axis, y-axis, 
z-axis, respectively. 


A relatively weak shear more strongly suppresses the concentration fluctuations along the shear gradi¬ 
ent direction. Thus, a perfect parallel phase is formed at the shear rate of y = 0.000005 in flgure[^(c). It 
is interesting that when we slightly increase the shear rate, the model transforms into a perpendicular 
lamellae with some defect as shown in figure [^(d). Then, continuously increasing the shear rate to an 
appropriate value, we get a concentric cylindrical structure where the block copolymer is surrounded by 
homopolymer in figure[^(e). Finally, at a higher shear rate, since the shear flow more strongly suppresses 
the concentration fluctuation along the vorticity direction, the model turns to the ordered perpendicular 
phase, as can be seen in figure |^(f). Besides this, the phenomenon also demonstrates that a low shear 
rate is beneficial for the formation of a parallel lamellar structure; and the lamellar phase preferentially 
adopts a perpendicular orientation at a high shear rate. 

As concerns the phenomenon under a simple steady shear flow, we can explain that it is a result of 
a competition between the shear flow and the concentration fluctuation. The imposition of a shear flow 
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suppresses the concentration fluctuation inordinately. Under a relatively weaker shear, the fluctuation 
effect is larger than the field effect, hence the chain cannot completely stretch along the shear direction. 
However, we have to admit that a weak shear also suppresses the concentration fluctuation to some 
extent, and it suppresses the fluctuation more strongly along the velocity gradient axis, i.e., y axis. This is 
the main cause of preferential adoption of the parallel orientation under a weak shear flow. The higher 
the shear rate is, the greater the field effect is. Under a stronger shear flow, the field effect is predominant, 
and it intensively suppresses the concentration fluctuation. Moreover, a stronger shear more strongly 
suppresses the fluctuation along the vorticity direction, i.e., z axis, thus, it transforms into the ordered 
perpendicular lamellar structure when the shear rate increases to a certain degree. This is in agreement 
with the previous study of the block copolymer under shear flow by K. A. Koppi et al. in experiments, (tT) 
and by Igor Rychkov in theoretical studies, (23, respectively. 

To consider this case in greater depth, we 
would numerically calculate the domain sizes 
Riit] {i = X, y,orz) in the x, y, or z direction as 
a function of time. The domain sizes can be 
derived from the inverse of the first moment of 
the structure factor S(k, t] as follows; 





(3.1) 


where 




J dkfc/S(k, Ri j^ dkS(k, t] . (3.2) 


In fact, the structure factor S(k, t) is determined 
by the Fourier component of the spatial con¬ 
centration distribution I521I55I . It is defined by 
S(k, r) = (|(/)(k, f)|>^ (^ . Figureshows the time 
evolution of the microdomain sizes Ri{t) in the 
X, y, and z directions as a function of time in the 
double-logarithmic plots. The results are averaged 
over ten independent runs. 

From figure (a), we can see that the mi¬ 
crodomain size Rx in equilibrium gradually in¬ 
creases (from curve a to curve e) as the shear rate 
increases, and reaches a maximum value in curve 
e which corresponds to the shear rate f - 0.0006 
in figure [^(e). Confinuously increasing the shear 
rate, Rx in equilibrium remains constant (from 
curve e to /). This fact shows that the growth of 
the microdomain size takes place along the x-axis 
(i.e., the shear flow) with an increase of the shear 
rate. We also see that the difference between the 
curves e and / is that a relatively short time is 
needed to reach the equilibrium at a higher shear 
rate. Although the change of the microdomain size 
Ry in figure (b) and Rz in figure (c) from the 
curve a to curve / is more complex, it is consis¬ 
tent with the figure complefely. We clearly see 
a steep increase from curves b to c and a sharp 
decrease from curves c to d in figure [^(b); on the 
contrary, we see a sudden decrease from curves b to c and a rapid increase from curves c to d in fig- 
ure[^(c). It means that we obtain the order parallel phase under the shear rate of curve c, and the defec¬ 
tive perpendicular phase under the shear rate of curve d. Moreover, ify(c) » Ry{f] and Rz(c] « Rzif) 


t 


Figure 3. Time evolution of the characteristic sizes 
of microdomains Rxit], Ry{t) and Rztt) for differ¬ 
ent shear rates, j - 0.0,0.000001,0.000005,0.000006, 
0.0006, and 0.001 in curves a-f. 


23801-6 
















Microphase transitions of biock copoiymer/homopoiymer under shear flow 


confirm that the lamellae adopt a parallel orientation corresponding to figure (c) and perpendicular 
orientation corresponding to figure respectively. From figure we also see that the microdomain 
size grows slower at an early stage while it changes apparently in the middle stage, then it is stable at the 
later stage. With respect to the lower shear rate, the higher shear makes the microdomain grow faster, 
untn achieving a constant value. 




1x10 


-1x10 

-2x10" 


— *—N1 

— N2 


0 200 400 600 800 1000 

shear strain 


Figure 4. The first and the second normal stress differences JVi and N 2 as a function of the shear strain 
for (a) 7 = 0.000005; (b) 7 = 0.001. 



Figure 5. The shear viscosity 77 as a function of the shear rate 7 for the block copolymer and homopolymer. 

Owing to the dual properties including elastic deformation and viscous flow, the polymer liquid often 
shows more complex properties under the external force. We know that it is meaningless to study the 
absolute value of the normal stress component individually, and the difference value of the normal com¬ 
ponent always stays the same, so we can use the first and the second normal stress differences value to 
describe the stress state. For a simple shear flow, we use the shear viscosity to describe the viscosity of a 
system. 

Figure]^ shows the first and the second normal stress differences Ni and N 2 as a function of shear 
strain under weak and strong shear, respectively. We can clearly see that Ni and N 2 are sensitively depen¬ 
dent on the direction of the lamellae. Figure]^ (a) corresponds to the case of forming a parallel lamellar 
structure. It indicates that Ni > 0, and it increases apparently along with an increase of the shear strain 
untn achieving the constant value. On the contrary, N 2 < 0, and it decreases apparently along with an 
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increase of the shear strain until reaching a constant value. Figure[^(b) corresponds to the case of form¬ 
ing perpendicular lamellar structure. Compared with figure[^(a), it indicates that ATi > A /2 > 0, and they 
increase along with an increase of the shear strain until reaching a constant value. 

As we know, the original formulas of the first and the second normal stress differences are Ni - 
Oxx - Oyy, and A /2 = Oyy - ozz- xx< Gzz ^fe the notmal principal stresses, and express the tension 
in X, y, and z directions, respectively. Here, we should consider the pressure that is imposed on the 
system in x, y, and z directions in order to correspond to the formulas <2.15) and <2.16) . Therefore, we 
have a good explanation that N\ and N 2 are sensitively dependent on the direction of the lamellae. We 
can take no account of the stress of the x-axis direction since x-axis is the direction of the shear flow, 
and it is mainly affected by the tangential stress. For the case of forming a parallel alignment, A /2 < 0, it 
means that the pressure in y-axis direction is smaller than that in z-axis direction; for the case of forming 
a perpendicular alignment, A /2 > 0 , it means that the pressure in y-axis direction is larger than that in 
z-axis direction. This verifies that the low shear rate is beneficial fo the formation of a parallel lamellar 
structure; and the lamellar phase preferentially adopts the perpendicular orientation at a high shear rate. 

Figure shows the calculated shear viscosity of the block copolymer/homopolymer as a function of 
the shear rate. Generally speaking, the shear viscosity gradually decreases with an increase of the shear 
rate, corresponding to shear thinning, in which the layer width is slightly decreased so that the transverse 
lamellae can adopt parallel or perpendicular orientation (35). However, in addition to the shear thinning, 
there also exists a shear thickening in figurej^ corresponding to the local convexes. The first larger local 
convex corresponds to the figure [^(b), the width of the domain is larger than the others; and the second 
local convex corresponds to the concentric cylindrical structure in figure|^(e) for the same reason as the 
first one. Points A and B in figurej^are the viscosity under the shear rates 0.000005 and 0.001, respectively. 
This indicates that the viscosity of the perpendicular phase is by far smaller than the parallel phase. 



Figure 6. (Color online) Simulation for block copolymer/homopolymer at different times corresponding 
to the weak and strong shear of 7 = 0.000005 and 7 = 0.001. The initial state is disordered, (ai): t - 50000, 

(hi): f = 500000, (ci): f = 5000000. Other details are the same as in figure 2. 

Moreover, in order to monitor the processes of forming parallel and perpendicular lamellar structure 
in detail, we show the simulation snapshots at different times for 7 = 0.000005 and 7 = 0.001 in figurej^ 
In other words, the two different shear rates represent the weak and the strong shear mentioned before. 

We simulate the occurrence of the parallel lamellar structure under weak shear flow in figures[^(al)- 
(cl). As could be expected, it is an irregular structure in the initial period of time, such as the structure in 
figurel^fal). With the time evolution, the domain morphology begins to tend to the parallel lamellae in 
figure 6 (bl). Then, the domain morphology gradually evolves into the parallel phase with some defects. 
Due to the shear flow, the defects become fewer and an ordered parallel lamellar structure appears in fig- 
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time steps 

Figure 7. The order parameter S of block copolymer/hopolymer as a function of time for weak and strong 
shear, corresponds to figure 6 . a: 7 = 0.000005 and b: 7 = 0.001. 


ure[^(cl). For comparison, we also consider the occurrence of a perpendicular lamellar structure under 
a strong shear flow. As shown in figures [^(a2)-(c2), the formation of lamellae phase only needs a very 
short time. At f = 50000, the system has formed a perpendicular lamellar structure with more defects. 
Moreover, compared with the disordered structure in figure[^(bl), it has formed a perfect perpendicular 
phase in figure [^(b2) because the shear field is so strong that it has an overwhelming superiority. Com¬ 
paratively speaking, the process of forming an ordered parallel phase is very slow because the external 
field is weaker and it has no advantage in competition with the concentration fluctuation. 

MeanwhUe, to verify the result further, we calculate the orientational order parameter S of the block 
copolymer/homopolymer under a weak and strong shear flow, 

S= <2cos^0-l>, (3.3) 

here, we choose to calculate the order parameter of yz plane. Thus, 6 is the angle between the unit 
normal vector of the lamellae in yz plane and the unit normal vector in the velocity gradient axis I72II73I . 
S is the order parameter, S = 0 represents the completely disordered state; S = 1 means a completely 
orientational phase parallel to the velocity gradient direction; and S = -1 means another completely 
oriented phase perpendicular to the velocity gradient direction. Figure[^displays the order parameter at 
two typical shear rates, 7 = 0.000005 and 7 = 0.001. 

As depicted in figure the initial value of the order parameter S in two cases is around zero. This 
means that the original state of the system is completely irregular. For the case of 7 = 0.000005 in curve 
a, the order parameter is in a mess before the time steps t - 750000, so the microphase structures are 
disordered in figures [^(al)-(bl). With time evolution, the order parameter S increases slowly because 
the shear flow gradually plays a certain role. Finally it reaches a certain value which suggests that the 
system forms an ordered structure, i.e., a parallel lamellar structure, which corresponds to figure[^(cl). 
In comparison with the above result, the order parameter S in curve b decreases quickly with an increase 
of time under strong shear, i.e., 7 = 0.001, and it drops the lowest in a remarkably short time. It corre¬ 
sponds to figures [^(a 2 )-(c 2 ), where the morphology quickly reaches the final equilibrium structure. In 
addition, the order parameter S has almost no change at the late stage, so the microphase structures that 
we finally obtain are stable. 

Next, we impose these two typical shear rates on the equilibrium system which arise from the self- 
assembly of block copolymer/homopolymer under zero-shear and explore the simulation snapshots at 
different times in figure [ 8 ] The equilibrium phase under zero-shear indicates that the initial state is or¬ 
dered as shown in figure]^ (a). We see that the domain structure merely tilts along with the shear flow 
under a weak shear flow as shown in figures [^(bl)-(cl). At the same time, a strong shear flow makes 
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Figure 8. (Color online) Simulation for block copolymer/homopolymer at different times under a weak 
and strong shear of 7 = 0.000005 and 7 = 0.001. The initial state is ordered, (a): f = 0, (hi): t = 1500000, 
(cl): f = 3500000, (b2): t = 50000, (c2): t = 3500000. Other details are the same as in figure 2. 


the ordered lamellar structure break up quickly and reform the perpendicular lamellar structure with 
some defects in figure |^(b2). The defects are eliminated with the time evolution, and it transformed to 
the perpendicular phase with some perforated lamellae in figure[^(c 2 ). 
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Figure 9. The order parameter S of block copolymer/hopolymer as a function of time for weak and strong 
shear, corresponds to figure 6 . a: 7 = 0.000005 and b: 7 = 0.001. 

Figureshows the order parameter corresponding to figure]^ We choose to calculate the order pa¬ 
rameter of xy plane. Thus, 6 is the angle between the unit normal vector of the lamellae in xy plane and 
the unit normal vector in the shear flow axis. S-0 represents a completely disordered state, S = 1 means 
a completely orientational phase parallel to the shear flow axis. In the case of 7 = 0.000005, S quickly de¬ 
creases to reach its minimum, and then slightly increases around S = 0. This indicates that although the 
shear flow is too weak to break up the ordered transverse, the weak shear could make the ordered phase 
turn to a disordered phase tilting to the x-axis. In the case of 7 = 0.001, S also quickly decreases to reach its 
minimum around 0 , but then it increases to a maximum around 1 in a short time and keeps a stable state 
during the next time. This result is consistent with the structural evolution shown in figures|^(a)-(c2). 
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From the above result, we can give the following explanations. When the initial state is ordered, 
the weak shear that we imposed is too weak to break up the ordered transverse phase. Thus, it cannot 
transform to another ordered structure. However, as figures [^(al)-(cl) mentioned, although the shear 
flow is also weaker, the initial structure is disordered, it is easy to form another ordered structure with 
time evolution. In addition, compared with the transformation from a transverse to a parallel phase 
of the block copolymer/nanorod composites (SHI. the whole interfacial tension increases when doped 
with the homopolymer, so the interfacial tension of the block copolymer/homopolymer is larger than the 
block copolymer/nanorod, and then the former’s order-disorder transition temperature (Todt) is higher 
than the latter’s. Thereby, it makes the order-order transition more difficult under a weak shear. In other 
words, this is the effect of the homopolymer. However, under a strong shear, the domain morphologies all 
finally turn to the perpendicular phase in the above mentioned cases including block copolymer/nanorod 
composites. This is because the shear flow is strong enough and absolutely superior, notwithstanding 
whether the initial state is ordered or disordered, and whether the system is block copolymer/nanorod 
or block copolymer/homopolymer. 

From another perspective, the thermal history cannot be ignored owing to the initial structure is the 
equilibrium phase under zero-shear. In fact, this is equivalent to the fact that the imposed shear flow after 
the phase is separated completely. So, we consider that thermal history also has a certain effect on the 
difference between figurej^and figurej^ Seung Su Kim and M. L. Cerrada (73 have studied the effect 
of thermal history on the phase behavior of polymeric materials, and it showed that some contradictory 
results in the references are likely to be caused by the different thermal history of a sample. 


4. Conclusions 

In this paper, the cell dynamics simulation is used to investigate the phase behavior and rheological 
properties of a block copolymer and homopolymer mixture subjected to a steady shear flow. Different 
morphologies corresponding to different shear rates, and the phase transitions occur from the transverse 
to parallel and then to perpendicular lamellar structure with an increase of the shear rate. It is a result of 
the competition the between shear flow and the concentration fluctuation. The weak shear more strongly 
suppresses the concentration fluctuation along the velocity gradient axis, while a strong shear strongly 
suppresses it along the vorticity axis, which is caused by the formation of a parallel and perpendicular 
phase under a weak and strong shear, respectively. This is in agreement with the previous work on block 
copolymer which was studied by K. A. Koppi et al. in experiments (71), and by Rychkov Igor in theory 
(2^ . Rheological properties being considered show that the first and second normal stress differences 
are sensitively dependent on the direction of the lamellae. The shear viscosity gradually decreases with 
an increase of the shear rate but there is a shear thickening at a certain stage. 

Moreover, we specifically explore the self-assembly of the block copolymer/homopolymer under a 
weak and strong shear in two different initial states, respectively, and obtain different results. The order 
parameter is in good agreement with this result. Interestingly, when the initial state is an ordered struc¬ 
ture, our results are just in agreement with the conclusion under strong shear studied by L. L. He 
Under a weak shear, the domain structure merely tilts along with the shear flow. This is because the whole 
interfacial tension of the block copolymer/homopolymer is larger than the block copolymer/nanorod, and 
the order-order transition is more difficult under a weak shear. In other words, this is the effect of the 
homopolymer. Certainly, with regard to different results in different initial states, the thermal history 
also has a certain effect. Therefore, it is indispensable to apply the shear field at an appropriate time if 
we want to get what we want. Our result could provide a guideline for an experimentalist, and the model 
system can also give a simple way to realize an orientational order transition in soft materials through 
imposing a shear flow. 
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MiKpo(|)a3Hi nepexoAi/i 6 /iomhi/iii Kono/iiMep/roMono/iiMep sa 

yMOB 3CyB0B0Y TeMM 

fO. fyo, /lx. llleHb, B. BeHb, X. By, M. CyH, /lx. fleH 

lilKO/ia xiMii i Marepia/iosHaBCTBa, neflarori^HUM yHiBepcuTer LLIaHbci,/liHb(|)eHb, 041004, Ki/iraw 

3 MOTOK) flOC/iiflxeHHfl cfiaBOBOi noBefliHKi/i cyMimi 6;iOHHnfi Kono;iiMep/roMono;iiMep nifl fliero cra/ioT 3cyBO- 
BOi TeHM BUKopucraHO MOfle/iiOBaHHfl fli/iHaMiKM KOMipoK. (PaaoBi nepexoflM BiflOyBaroibCB Bifl nonepenHOi flo 
napayie/ibHOi, a noriM flo nepneHfli/iKy/iflpHOi /laMoyiapHoT crpyKTypi/i 3 pocTOM acyBOBoT ujBUflKOCTi, Lpo e pe- 
ay/ibiaroM B3aeMOfliT Mi>K acyBOBOK) Tenieio i cfi/iyKTyapiBMn KOHpeHTpapiT. Vci peo;iori4Hi B/iacTHBOcri, raxi ak 
piaHi/ipi HopMa/ibHoro HanpyxeHHa i acyaoBa B'^axicTb, ricHO noB'aaaHi a HanpaMKOM /laMe/in. Bi/ibiue roro, 
aoKpoMa HaMU flOc;iiflxeHO cfjaaoBy noBefliuKy i napaMerp nopBflKy npn cyiaOxoMy i ci/i;ibHOMy acyeoBi flBOx pi- 
3HI/IX noHaiKOBnx craHiB, a BpaxyBaHH^M npn pbOMy Ba>iaii/iBOCTi TepMiHHOi icropii. fl/ia OAep>KaHHfl OaxaHoro 
peay/ibiary acyBoae no;ie c/iifl npuK/iac™ y BiflnoBiflHnii momomt 4acy. Peay/ibiarn piei po6oTi/i aaOeanenyiOTb 
yierKMM cnoci6 CTBopeHHB BnopBflKOBaHnx OeafleefieKTHi/ix Marepia/iia bk b excnepuMeHTi, rax i b iHxeHepHm 
TexHO/ioriT aa paxynox BBefleHHfl acyBOBOiTeHu. 

KnioMOBi cjioBa: caMOCKyn^eHHfi, 6no^Hi/i£i KononiMep, roMononiMep, acyeoBa re^ifi 
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